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Abstract 



One-loop corrections to the fermion rest mass Mi , wave function renormaliza- 
tion Z2 and speed of light renormalization Co are presented for lattice actions 
that combine improved glue with clover or D234 quark actions and keep the 
temporal and spatial lattice spacings, at and as, distinct. We explore a range 
of values for the anisotropy parameter x = o-s/o-t and treat both massive and 
massless fermions. 
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I. INTRODUCTION 



Examples of successful employment of anisotropic lattices in lattice QCD simulations 
have been increasing lately. They include extensive studies of the glueball spectrum [|l|, 
investigations of heavy hybrid states and calculations of quarkonium fine structure . 
In most cases one is dealing with large states requiring large spatial volumes and also signals 
that can only be extracted from high statistics data. Working with highly improved actions 
on coarse lattices helps with the large volume and statistics problems, however, a coarse 
temporal lattice spacing means that correlation functions fall off very rapidly. This last 
problem can be circumvented by going to an anisotropic lattice which allows for a much 
finer temporal grid. The correlation functions can be sampled much more frequently in a 
given physical time region where the signal is still good. 

Another potential use of anisotropic lattices would be in simulations of matrix elements 
in hadronic states with large momenta. These typically occur in semileptonic decays of 
heavy hadrons. Once one goes beyond spectrum calculations to matrix elements, one is 
faced with the matching problem between operators in the coarse highly improved lattice 
theory and continuum QCD. In this article we take a first step in accumulating necessary 
renormalization factors based on perturbation theory. We carry out the one-loop self energy 
calculation in several anisotropic improved quark actions. This gives us renormalization of 
the rest mass Mi, the wave function renormalization and the "speed of light" renor- 
malization, Co, a quantity which will be defined more precisely in the following sections. 
We treat both massive and massless quarks. Perturbation theory can be used not only in 
operator matchings, but also to fix parameters in the lattice actions. Co is an example of 
one such parameter. 

In the next section we introduce the gauge and quark actions considered in this article. 
Section 3 describes the general formalism that we employ for the self energy calculation 
with massive fermions. We follow closely the work of the Fermilab group |^ which we 
could straightforwardly extend to anisotropic actions. Section 4 discusses specific one-loop 
contributions for mass, wavefunction and speed of light renormalizations. Our results are 
tabulated in section 5 for various choices of actions, fermion masses and degree of anisotropy. 
Some calculational details are left for appendices, where we describe Feynman rules and IR 
subtractions in our calculations. 



II. GAUGE AND QUARK ACTIONS 



We work with two classes of gauge actions denoted Sq and S]^ with 
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The X sum is over lattice sites and the variable s runs over spatial directions. (3 = 2Nc/g'^, 
X is the anisotropy parameter 



and 



P, 



X = as /at 



Real (Tr{Uf^{x)U,y{x + af^)Uj^{x + a^)Ul{x)} 



(3) 



(4) 



Rtiu = -Tj-Real (Tr{U^{x)U^{x + a^)Uy{x + 2a^)[/^(x + a^ + ay)Ul{x + a^)Ul{x)yj . (5) 

Us and Ut are the tadpole improvement parameters Uq for spatial and temporal link variables 
respectively 

The parameters Cq and cf in action Sq are constrained to satisfy Cq + 8cf = 1. The 
Symanzik improved gauge action, in which 0(a2) errors are removed, corresponds to Cq = 
5/3 and cf = —1/12 0, whereas Cq = 3.648 and cf = —0.331, for x = 1; leads to one of the 
RG improved Iwasaki actions [0]. In Sq parameters have been fixed to the Symanzik values. 
We will be working mainly with Symanzik improved actions and present RG improved results 
only for a few cases. We note that the action Sq is intrinsically asymmetric even for the 
isotropic limit x = 1- 

The most highly improved quark action that we have analysed is the D23A action 

.3, 
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The quark fields \l/c and the dimensionless lattice fields '^l are related through 
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The dimensionless derivatives V^"-* and field strength tensors F'^'^ are tadpole improved |§] 
and defined in the Appendix. We use the convention cr^j, = ^[7^,7^] and set r = 1 in all 
our calculations. At tree-level the coefficients Cq, C3, Cst, C4, and Cp are equal to 
one. The quark action is then tree- level accurate through 0{a?^) and O(a^). Cq is what 
we call the "speed of light". This parameter is adjusted, in general either perturbatively 
or nonperturbatively, to ensure correct dispersion relations for particles. In anticipation of 
working on anisotropic lattices with at much finer than a^, one can drop the higher order 
improvement terms in the temporal derivatives by setting Cst = Cu = 0, without loosing 
accuracy. We call this action iS|,234- 



E l^tVt + — 7 ■ (V - ICsV^'^) + atmo 
X [ X 6 
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The familiar 0{a) accurate clover quark |11| action corresponds to setting C3 = C4 = in 
the above and using a less improved field strength tensor F'^'^ (also defined in the Appendix) 
rather than F'^'^. 



Sdover = J2^lI It^t + —7 ■ V + atmo 

X \ X 
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We have carried out one-loop self energy calculations for several combinations of the 
above gauge and quark actions, for both massless and massive quarks. We list the specific 
actions considered in Table I. For actions and massless results have already appeared 
in [|12|. We agree with their results and we include these cases here for completeness. With 
action S'" we treat only the massless case, since our formalism for massive quarks, following 
0, requires that the only time derivatives be in the Vt and v[^^ terms. 
Sdover satisfy this condition, but S{)2m does not. 



Both iSl/234 and 



III. GENERAL FORMALISM FOR SELF ENERGY CALCULATIONS 

In this section we summarize the formalism for self energy calculations, along the lines 
of reference 0. Perturbative calculations for massive Wilson quarks are also described in 
reference ||TB[ . We concentrate on the massive case, since massless lattice perturbation theory 



has been in the literature for decades. 

For massive fermions we use either iS|,234 Sdover- The fermion self energy S(p) is 
defined in terms of the momentum space propagators G{p) and Go{p) for the full and free 
theories respectively, as 
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g-\p) = g-\p)-j:{p). (11) 

Carrying out the Fourier transform in po one defines 

^ Z2ip)e-^^P^%ro, + ... . (12) 

Tproj is a projection operator in Dirac space. The elhpses refer to lattice artifacts and 
additional multi-particle states that could be created by the lattice fermion field operator \E' 
beyond the single quark state. The rest mass, Mi, is defined as 

Ml = E{p = 0). (13) 

We do not consider the kinetic mass, M2 in this article. We will renormalize at the point 
{PojP) = {'iMijO) and define the wave function renormalization constant 

Z2 = Z2{p = 0)- (14) 
For a zero spatial momentum quark propagating forward in time one expects (t > 0) 

G(t,0)= ^e^Po*^(po,0) 

J-n/at ^TC 

= Z2e-^^^*i±^ + . . . . (15) 

Our goal in this section is to relate Z2 and Mi to parameters in the action and to S(p). In 
order to orient ourselves, however, it is useful to first consider the free case with S(p) = 0. 



A. Free Anisotropic Propagator 

The free propagator Go{po,p = 0) for both actions iSl/234 and Sdover becomes (for r = 1) 



— Gq{pq,p = 0) = — — ^— 

at z7o sm(atpo) + 0-1^0 + x - X cos(atpo) 

-i7o sin(atpo) + ajmo + X - X cos(atPo) 



(16) 



(sin(atpo))^ + [{atmo + x) - X cos(ajpo)]^ ' 
In terms of the variable 

z = e*"*f° = e"'^'-^ (17) 
{po = iE), one finds two zeros of the denominator corresponding to positive energy solutions. 



{atmo + x) - V {atmo + x)^ + 1 - X^ 
= (18) 
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and 



{atniQ + x) - J{atmo + x)^ + 1 - ^ ^ 

^' = ■ 

The other two zeros, Z2 and Z2 correspond to negative energy solutions, Z2 = 1/zi and 
Z2 = l/zi- The integral over pq in (p!5D can be done as a contour integral around the unit 
circle in the variable z. One picks up contributions from both positive energy solutions (for 
t > 0). 

/ , (I + 70) e-<^ 
pole at Zi : , (20) 

2 ^(atmo + x)^ + 1 - 

po/e at 2:1 : ^=, [21) 

2 y^(aimo + x)^ + 1 - 

with 

aiMj°) = - ln(^0, aM'^ = - In(zi). (22) 

Clearly, is the physical positive energy solution. The second solution zi is a lattice artifact, 
similar to the time doubler for r 7^ 1 in isotropic actions. The solution zi disappears in the 
isotropic limit, x I5 where ajM^^^ — ^> 00. In the same limit the physical solution zi goes 
over into the well known result 

zi i . (23) 

1 + atiriQ 

The gap between Mi^^ and m[^\ measured in units of l/a^ is 

a,(Mf) -Mf)) = x In (24) 

independent of mo. This becomes 00 at x = 1 and approaches 2 as x 00. The size 
of this gap, agSE > 2, is hence equal to or larger than the amount by which conventional 
spatial doublers are raised through the Wilson mechanism. We will henceforth ignore zi and 
concentrate on the physical pole at 2; = 2:1. Comparing (^) with (|15|) one sees that there is 
nontrivial mass dependent wave function renormalization even at tree-level with 

^(0) = I = I (25) 

' '[atmo + xy + l-X^ X sinh(atMf )) + cosh(atMf ^) ' 

This has been pointed out several times in the literature [0,^. 
A useful way to rewrite (|18D is 

a^mo + X = X cosh(atMf ^) + sinh (a^Mf^). (26) 
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B. Mass Renormalization 



In the interacting case one has a nontrivial which we write as 

atT.{p) = ijoBo{p, mo) sin(aiPo) + i- Y^lljBjip, mo) sm{asPj)] + C{p, mo). (27) 



The p = propagator becomes 



— G(po p-0) - ~^'")'o(^ ~ -^o) sin(afPo) + Ofmp + x - C - x cos(atPo) 
at °' (1 - Bo)2[sin(atpo)]2 + [(atmo + x) - - xcos(atpo)] 



where, Bq = Bq{j>q, ttiq) and C = C{po, mo) are evaluated at p = 0. If po = iE is the location 
of a pole in (|28|) , the following implicit equation must be satisfied. 

(1 - Bo{iE,mo)) sinh(ai_E') = ±[(atmo + x) - C{iE,mo) - xcosh{atE)]. (29) 

One can check that the "+" sign leads to the pole m}°^ in the free limit. Hence, the implicit 
equation for Mi is given by 

X cosh(afMi) + sinh(aiMi) = atmo + x + Bo{iMi, mg) sinh(atMi) — C{iMi,mo). (30) 

In a perturbative calculation of Mi one expands 

Ml = Mf ^ + a,Mf ^ + 0{al). (31) 

Bq and C in (|30|) start out 0(a;s), so through one- loop their argument can be replaced by 
the tree-level Mi^\ Expanding the LHS also through 0{as) and taking (p6[) into account, 
one finds 



(1) _ 5o(zMf\mo)sinh(aiMj°Vc(2Mf\mo) 



Xsmh{atM{') + cosh(atM| 
'(0) /(7o + l) _,.^{o) 



= _ tr ^^^L_^atJ:{po = iMr\p= 0)| , (32) 

where the trace is taken over Dirac space. We note that in the Mi'''* = 0,mo = limit, the 
7o part of the trace tr{(7o + 1)S} does not contribute and one has 

a,atMf^(0) = - tr{ai S(0)}/4 = -C(0,0). (33) 

In order to have massless quarks remain massless under renormalization, one needs to carry 
out additive mass renormalization and m'^^ in ( P^ ) requires a subtraction. This subtraction 
must be done without jeopardizing the pole condition (|30|) . There is a standard way to 
accomplish this. Let mc be the value of the bare quark mass parameter mo for which the 
physical quark rest mass vanishes (Mi = 0). Eqn.(|30|) then tells us that rric is implicitly 
defined through 



7 



atTJic - C(0, rric) = 0. (34) 

In equations such as p^) or (|30D one always has the combination a^mo — C. Using (^4|) one 
can add and subtract atrric so that 

atrriQ — C ^ dtirriQ — nic) — (C — C(0, rric)) = atm — C . (35) 

Previous derivations go through with mo replaced by 

m = niQ — rric (36) 

and C{iMi,mQ) by 

^(iMi, mo) = CiiMi, mo) - C(0, m^). (37) 

In most lattice simulations, iric and hence also m are determined nonperturbatively from 
the simulations themselves. For C, however, one still often uses the one-loop result 

^(iMf \ mo) = C{iM[^\ m) = C{iMi°\m) - C(0, 0). (38) 

is now given in terms of m rather than mo- We will be presenting our results as 
functions of asMi'\ with the understanding that the shift mo — > mo — rric has been carried 
out and that, for instance, Mi'^ is given by 

a^m + x = X cosh(atMf ^) + sinh(atMf ^) (39) 



rather than by p^). In (^) one needs to replace C by C. Our final formula for the one- loop 
mass correction, measured in units of l/a^ then becomes 

(1) _ So(^Mf \ m) sinh(atMf )) - ^(^Mf \ m) 

O^sO-s-ly^i.ub — X ~ 



Xsinh(atMf^) + cosh(atMf^' 



) tr 



(70 + 1) 



a,S(po = «Mf ,p = 0,m)-a,S(0,0,0) \. (40) 



This expression vanishes automatically for m}°^ = m = 0. We prefer to measure dimension- 
ful quantities in terms of l/a^ rather than 1/at. When exploring x > 1 it makes more sense 
to fix as and let at be arbitrarily fine, rather than to fix at and let a<j become arbitrarily 
coarse. In the isotropic limit (^) agrees with formulas in the literature [§,0. 



C. Wave Function Renormalization 

In order to extract a general formula for the wave function renormalization Z2 we need 
to find the residue of G{po,p = 0) at the pole po = iMi. In terms of the variable z the 
Fourier transform in (p!5| ) has the form 

(41) 



8 



where the integral is taken over the unit circle. To find the residue we expand the denomi- 
nator around Zi = e""'*^^^ 

\ / Z = Zl 

The contribution from the physical pole to G{t, 0) is then 



\ J \ J / Z = Z\ 

One finds for the numerator 

g{z = zi) = (70 + 1)(1 - BoitMi, m)) sinh(atMi) (44) 
and for the denominator 
2(1- Bo{iMi,m)) sinh(atMi) x 
Xsinh(a4Mi) + cosh(atMi) + ( ^^^~j[^^o(Po, m) sin(atPo) + C{po, m)] ) } (45) 



PQ=iM\ 



usmg 



urn --^{^^ . (46, 

One can now read off Z2 and after recognizing the last term in (^51) as derivatives acting on 
different parts of atT,{pQ,p = 0,m), one obtains 

Z^^ = xsmh{atM,) + cosh{atMi)+itr ( j ^E{po,p = 0,m) ) . (47) 

V 4 dpo 

The one-loop approximation to Z2 is obtained by expanding Mi once again in ctg. 

^2"^ = Xsinh(atMf ^) + cosh(atMf ^) + a.atMfJ^^ (x cosh(atMj°^) + sinh(atMj°^)) 



1 + ^a,M« (xcosh(a,Mf )) + sinh(a,Mf ))) 



4 dpo" 



+^trf^^^As(p„,p-'=0,m)) Zf] + 0(a,2). (48) 



In the last expression we have found it convenient to factor out the tree-level z'2^ ^ . Equa- 
tions (|47D and ( ^8| ) go over into the formulas of |^ in the isotropic limit. 
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D. Speed of Light Renormalization 



In order to discuss renormalization of the speed of light one needs to look at the inverse 
momentum space propagator at small but nonzero spatial momentum. 

O'tG = atGQ - atT.{p) 

= i'jo{l - Bq) sin(atpo) + i - Ylbi (^o^j - Bj) sm{asPj)] 

^ j 

+atm + X - X cos(atpo) - C, (49) 
with Kj = 1 for Sdover and Kj = (4 — cos(asPj))/3 for One can rewrite G ^{p) as 

atG~^{p) = (1 - Bq) J i7o sin(atpo) + i - Ylbi^^l^^—^r^ ?>in{asPj)] 

\ X j ^0-- Bq) 

+atm + X - X cos{atpo) - G. (50) 

Gq is adjusted so that for small asPj the relative coefficient of the 70 sin(atpo) and the 
•yjasPj/x terms remains equal to unity. {G^Kj — Bj) sin(asPj) = (Co — Bj)asPj for all quark 
actions in this limit (of course, Kj sm{asPj) is a better approximation to the continuum agPj 
in the D23A action than in the clover action), and one has 

Go = l + Bj{Go) - BoiGo) ^ 1 + B,{Go = 1) - Bo{Go = 1) + 0(«^). (51) 

Just as with Z2 we will define the speed of light renormalization at the zero spatial momen- 
tum mass shell point p = (iMi,0). From ( P7| ) the two terms Bj and Bq needed for Gq at 
one-loop can be extracted through 



4 \" diasPj] 



l tr(7oafS(^Mf\o)) 

Bo = -7 M m > (53) 

4 sinh (a^Mf)) 



or 

-7^n70 7^S(p)) m = 0. (54) 

/ p={0,0) 



4 V 9Po' 



We note that in the massive case there is nontrivial renormalization of Gq even in the isotropic 
limit X = 1) due to our noncovariant mass shell condition, p = {iMi,0). Nevertheless, we 
believe the above definition of the renormalization of the speed of light is a sensible and 
physical one. 



IV. ONE-LOOP CONTRIBUTIONS TO S(P) 

In the previous section one- loop corrections for Mi, Z2 and Gq were determined in terms 
of traces over S(p) or over derivatives acting on E(p). In this section we describe the lattice 



10 



perturbation theory diagrams that contribute to at one loop. For all quark actions 
considered one can write 

S (p) = TT's (p) + E*"'^ (p) + S*-^- (p) . (55) 

E'^^^ is the regular "rainbow" diagram, the only diagram that exists in a continuum self 
energy calculation. E*""^ denotes contributions from the lattice artifact tadpole diagram and 
j^t.t. QQj-Qgg fi-om perturbatively expanding the m^'s and Uts entering definitions of tadpole 
improved derivatives (see Appendix). The main idea behind tadpole-improved perturbation 
theory p is to have S*'*' cancel the bulk of E*"^*^. We find in many instances, especially with 
^clover, that cancellation is complete if one uses the Landau link definition for Ug and Ut and 
works in Landau gauge. 



A. S''^^?(p) 

In terms of the gauge propagator D^j,, quark propagator Gq and the vertex functions 
V/j,, one can write atS'"^^(p) as the following loop integral over the dimensionless momentum 
variables — vr <k^<ix. 

= 9^^^J o-V - k) ^^^^ — — y^iap - k, ap)| D^^{k, Og) 

(56) 

where, ap stands for {atPo, agp) and A' sin for Kq sin((ap — A;)o) or for Kj sin((ap — k)j)/x- 
Kfj_{k), Q{k), V^{k', k) and Dfj_^{k, ag) are detailed in the Appendix. The argument ag in the 
gauge propagator comes from the gauge fixing term, with ag = 1 and ag = corresponding 
to Feynman and Landau gauges respectively. Our codes have been written for general ag 
and we have used gauge invariance of Mi and Co as one check on our results. 



Equation (|56D has the familiar form for a self energy integral. The only subtlety is to verify 
that one is indeed calculating S''^^ measured in units of l/a^, given the conventions in our 
Feynman rules. As explained in the Appendix, we choose to work with a dimensionless 
momentum space gauge propagator -D^;^. It comes from the Fourier transform of the di- 
mensionless correlator ((a^A^) {a,jA,j)). The relation between D/j^^, and a more conventional 
propagator of dimension 1/ (energy)^, denoted D^u, is 

D,. = ^D,^. (57) 

Our vertex functions, V^, are also subtlely different from those in isotropic lattice perturba- 
tion theory. They keep tract of the 1/a^ in the derivatives, i.e. of whether one has a l/a^ 
or 1/at there. If are vertex functions normalized such that — >■ —27^ for all /i in the 
continuum limit, then the relation to the V^i of (^) is given by 
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Using (|5^) and ( |58D one can write 

,4^ r d^k 



at ~ \ Go ( at ~ 



4 




at Vfliy 



(2,r)« 




(58) 



(59) 



To evaluate ( |5B| ) we made extensive use of the symbolic manipulation package Mathematica. 
The integrals themselves were done using the VEGAS program The various steps 

involving Mathematica were to 1. calculate the products V^[— ^7 ■ -fCsin 2. carry 
out the Dirac traces such as tr{(l + 7o)S} ; 3. take derivatives with respect to external 
momenta ; 4. put things on the mass shell p = (iMi, 0) ; and 5. use trigonometric identities 
to re-express the full integrands in (|56| ) in terms of powers of A;^ = 2 sin(A;^/2). The last step 
facilitated speedy evaluation of the integrand by VEGAS. 

Both Ml and Co are physical quantities. In addition to being gauge invariant they are also 
IR finite. The wave function renormalization Z2, on the other hand, is gauge dependent and 
also generally logarithmically IR divergent. In any calculation of a physical quantity this IR 
divergence will eventually be cancelled by vertex corrections and/or matching to continuum 
operators. In this article we will isolate the gauge dependent IR divergence in Z2, verify 
that it is the same as in the corresponding continuum theory and present results for the 
remaining IR finite parts. The IR divergence is found in the contribution from Y7^^ to Z^- 
More specifically it resides in the following term in (|48|) 



(60) 



We adopt the method of reference |T3| to subtract IR divergent contributions inside inte- 
grands and rewrite ( |60|) as 



(70+1) ^ 

4 c/po 
F(meif,A,A), 



V — V 

at 



(0) 



-D^jy Z2 



(0) 



with 



and 



F{meff, A, A) = ^ J^sub{k, TUeff, A, A) 



4 r d^k 



3 J (27r)4- 



(61) 



(62) 



(63) 
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Explicit forms for J-'suh{k,mejf, A, X) and F{meff, A, X) are given in the Appendix. A is a 
gluon mass introduced to regulate IR divengences. J^sub has been constructed to match the 
same IR divergence as the first term inside the integral in (0). As a result the integral 
becomes independent of A. The other condition on J^sub is that the integral (^) be easy 
to do analytically. The simplest approach is to use a continuum self energy expression for 
J-'sub with an appropriate choice for the mass parameter nieff. The need to adjust niejf to 
optimize matching of the IR behaviours in J^sub and the lattice integrand, was emphasized 



in reference |]T3[ and following that work we find 

atnieff = sinh(atMf V"""^'^"^'^^ \ ' ^-^"^^^"^^^^-^^ ' _ (-54^ 



^ cosh(atMf ^) + X sinh(atMf 
1 + xsinh(atMf^^ 



The same J^sub and mf,ff work for both the clover and D234 quark actions since the IR 
structure of the two theories agree. Finally, A < vr in the above expressions is a cutoff 
imposed on J^sub so that Tsub = for /c^ > A^. The full expression (|6T|) must be independent 
of A. 



B. E*'^'^(p) 

The second contribution to is the tadpole contribution coming from the 

two-gluon emission vertices listed in the Appendix. For quark action 5^^234 has 

atT!""^{p) = ^ [^70 sin(atpo) - X cos{atPo)] A 



''00 



(3 + C3) sm{asPj) - 2C3 sin{2asPj) cos(y ) 



(3 + C4) cos{asPj) — C4 cos{2asPj) cos^(-^) 



(65) 



S* (p) in the case of Sdover is obtained by setting C3 = C4 = in the above expression. 
The appropriate traces and derivatives with respect to external momenta can be carried out 
immediately and one has 



i[sinh(aiMf )) + X cosh{atMi'^)] J, Doo + ^ E, Ik D. 



(0)^ 



2x ■ 



"J3 



s. 



clover 



i[sinh(a,Mf )) + X cosh(a,M|"^)] /, Doo + ^ L %[4 - cos^^)] 5^^234 



r(0) 



(66) 



,,J(7o + l)As-'^^ 
[ 4 dpo 



p=(iMf\o) 



--[cosh(atMi^ 



(0)^ 



Xsinh(atMf))] / Dqo (67) 



II 



Sclover <fc '^D234 
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{2 -fk ^jj 2 -fk -^00 ^clover 

(68) 

where B*"'^ and Bq"-'^ are the contributions from the tadpole diagram to (|52|) and (^) or 
(^). All the integrals are IR finite and very easy to carry out numerically. Contributions 
from S*"*^ typically dominate over those from E*"*^^ but the bulk if not all of it is cancelled 
by E*-*-. 



The lattice covariant derivatives in the quark actions are tadpole-improved. They are 
listed in the Appendix. In momentum space one has, for instance 

i sin(a^p^)/M^ ^ i sin(a^p^) [1 + asW^f^] + O(a^), (69) 

where we have perturbatively expanded 

= + 0{al). (70) 



Even in the absence of the regular and tadpole one-loop diagrams there are hence 0{as) 
terms in the quark propagator. We denote the inverse quark propagator with the m^'s still 
in place as ^0,^0(^)5 so that Gq (p) = Gq^^q^^^p). Through 0{as) eqn. ([TT|) can be written 

as 

G'^ = Go;io - S""^ - S*'*'^ = Gq^- - ^""^ - S*-*- (71) 

or 

S*'*' = Gq — Gq Q. (72) 



From the difference in ( [72l) one sees that one link hops bring in factors of 1 — I/m^j ^ —agU 
and two link hops factors of 1 — 1 /m^ ~ —2asU^^^ etc. Using these rules one finds 

at (p) = [-z7osin(aiPo) + Xcos(atpo)] + 

(Xsuf^ — {-hj [(3 + C3) sm{asPj) - C3 sm{2asPj)] 
^ j 

Ui) cos{asPj) — 

The relevant traces and derivatives become 



(2) 



G 1 

+ [(3 + C4) cos{asPj) - cos{2asPj)] | . (73) 
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—tr 



(7o + l)„ w..." 



p={iM[°\6) 



-[sinh(atMf^) +xcosh(atM| 



(2) 3 



(2) 



5, 



(0)^ 



,(2) 



^ 2 "s^s 



clover 



all 



(74) 



1 4 dpo 



cosh(atMj°^) + X sinh(atM{''0] 



(0)^ 



,(2) 



5, 



clover 



II 

D234 



(75) 



5] 



i.i. 



- fnf )) 



The Landau mean link definition of m„ is given by 



(2) 



1 



S, 



clover 



qll 



'-ID 

2Jk 



0). 



(76) 



(77) 



If one evaluates S*'^'^ in Landau gauge then (|6D & (|7D , (l^) & (|7|) and (|68|) & (|7|) cancel 
for S clover- (§3) ^ (0) caucel for 5^/234 and for the other two traces cancellation is 
almost complete. The difference between contributions from S*""^ and S* *' would go away 
if one replaces cos^(A;/2) and sin^(A;/2) by their averages 1/2. Hence, it is easy to see in 
this calculation how tadpole improving terms in the lattice action eliminates lattice artifact 
contributions in perturbation theory. 



V. RESULTS 

In this section we summarize results for the one-loop coefficients, OsM^^j^j,, and C'^^ 
for mass, wave function and speed of light renormalizations respectively. These follow from 
equations (|40|) , (^81) and (0) - (|54[) and each has, as explained in the previous section. 



contributions from regular, tadpole and t.i. diagrams. The numbers in our Tables are 
coefficients of a^. The Landau mean link definition of is used throughout to implement 
tadpole improvement. 



In Table 11 we present results for a^Mj-^j^^ for action for several values of asMi'\ We 
list separately contributions from Y7^^ , E*"'^ and S* *'. The fourth column gives the gauge 
invariant combination (reg + tad) and the sixth column gives asM^^}^^^^^ = (reg + tad + 
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t.i.), the full tadpole improved one-loop correction before subtraction. Carrying out the 
subtraction according to (|40|), one obtains a^M^^^J^^ which is given in the last column 



OsJWi.suf, - — — — —j-r — — — T^T-. {(a) 

Xs mh(atMr) + cosh(afMr^ 



All our calculations have been carried out for two values of the gauge fixing parameter a^, 
1.0 and 0.0. Table II lists both sets of results and one sees that gauge invariant quantities are 
independent of Ug within numerical integration errors ( which we take to be at the ±0.003 
to ±0.006 level depending on the mass). Our results for agMi'^ = agree with those from 



We plot ttsM^^j^fj versus agMi'^ in Fig. 1 . One sees that the mass dependence is smooth and 

that one reaches saturation rapidly already around OsMj-"^ ~ 3.0 — 5.0. We also compare 
with non-tadpole improved results for which the curve saturates around 1.827 rather than 
around 1.077. In considering the large mass limit it is useful to note that the factor Z2^^ = 
l/[xsinh(atMf^) ± cosh(atMf^)] appearing in P0| ) and ( ^H]) vanishes exponentially in this 
limit. The only terms that survive into the static limit are those where Zg'^'* is multiplied 
by an exponentially increasing function of asM[^\ It is easy to see, for instance, that the 
subtraction term in ( ffHl) or the spatial tadpoles in (|6^) become irrelevant in the static limit. 
Furthermore OsM^^j^j, becomes identical for Swuson, <Sdover and Sj^234: this limit and the 
only difference between the current calculations and those of reference resides in the 
glue action (we have verified that by switching to the unimproved Wilson glue action results 
of [|l^ are reproduced). 

In Tables III, IV and V we summarize results for OsMj^^j^^ for the other actions listed in Table 

I. We also list the combination (reg ± tad) for each case. If one chooses to implement tadpole 
improvement differently from what we have done here or decides not to tadpole improve, then 
agM^^j^f^ can be calculated straightforwardly from (reg ± tad) and the formulas presented 
in this paper. For action S'" we list only massless results for reasons explained in section 

II. The dispersion relations of this and similar actions are discussed in reference |T^. Our 



choices for anisotropy values in Tables IV and V, were dictated in part with an eye towards 
practical numerical simulations. Values such as x = 3.6 and x = 5.3 were taken from recent 
work on nonperturbative determinations of the renormalized anisotropy in pure glue theory 
|TP|. For one value agM^^^ = 1 we plot agM^^j^^ versus x iii Fig- 2 using the action S^. 
From Tables IV and V and from Fig. 2 one sees that the dependence of a^M^^j^j, on x is 
very mild. Effects of tadpole improvement are significant only for small values of x- This is 



because due to cancellations in the subtraction of equation (^8]) only the temporal tadpole 
and the temporal Landau link term u[^^ contribute to ctsM^^^j^^ and both these become small 
as X increases. 

Finally we mention that the one-loop expression for the critical bare mass rric is given by 

amc = --a.Mg_,(0) ± Oial). (79) 
X 
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B. 



Starting with ( ^8]) we define 



Z2 = zf [1 + 



(1) 



+ 0{al)] 



with 



3^ [-2 + 



[Ola 



ln(A2) 
hi(A2) 



m > 



(80) 



(81) 



A is the gluon mass in units of l/a^- It is the coefficient of after factoring out that 
has the same IR ln(A) and In(am) structure as the continuum wave function renormalization 
constant. From (^) one also sees that there are two contributions to Z'^\ one coming from 
the d/dpQ derivative term and the second from the expansion of Mi. Accordingly we write 



7(1) 



, 7(1) 



(82) 



In the literature ^'g^L is not always included as part of the definition of Z. 



(1) 



4!io alone 

with unimproved Wilson glue goes over in the large mass limit to the wave function renor- 



malization of reference [|T^]. Including Z2^l^^ leads to the static result of reference |^ which 



has been used in many subsequent static calculations, for instance in This latter static 
value also corresponds to the large mass limit of the one-loop Z2 calculated in many versions 
of NRQCD actions H- 



Table VI presents results for -Z'g^jp^ and the full ^2^-* for the action S^. Again we agree with 



reference \12] for m = 0. However, one notices that the massive data do not tend towards 



the massless result as asM^'^^ decreases. This is because our massive numbers include In(am) 
contributions which will eventually diverge, whereas in the massless theory we have set the 
fermion mass identical to zero from the beginning. This leads to different IR structure for 
the two theories (see Appendix B for some further discussions). In a matching calculation 
one will be looking at differences between the lattice and continuum Z2. As long as IR 
divergences are handled in the same manner in the lattice and continuum evaluations, one 
should not run into any problems and the m — >■ limit should be smooth. For instance, 
using dimensional regularization in the MS scheme one finds in Feynman gauge the UV 
finite continuum results 



,7(1) cont. 
^2 



Stt 

J_ 

Stt 



m 







ln(^)+21n(f )-4 



(83) 



m > 



Taking the difference between continuum and lattice wave function renormalization con- 
stants, it makes sense to consider the following subtracted Z2 factors. 



7(1) 



7(1) 



J_ 1 

Stt 2 



m 



J_ 

Stt 



[31n(a,Mf^)2-4 





m > 



(84) 
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Numbers for Z\ j^jj- are given in Table VII and one sees that the m ^ behaviour is smooth. 

In Tables VIII through XII we present ■^2,ipo ^2 ^ other actions. One does not 
find any dramatic changes with differing actions and/or anisotropy. The IR subtractions of 
Appendix B worked well in all actions for a^M}"^ < ~ 5. For larger masses VEGAS errors 
became large especially for x = 1. Hence, we only present results up to asM^'^ = 5. For 
X > 1 problems were less severe in general. In Tables VI - XII the numerical integration 
errors are at the ±0.02 level for agMi*^ = 5.0 and x < 3, at the ±0.006 level for asMi^^ = 
0.01, 0.05 and 0.10 and at the ±0.004 level for all other cases. A more sophisticated method 
for handling IR divergent integrals appears necessary if accurate results are required for 
larger masses. Many quantities, however, are close to saturation by the time one reaches 
a,Mf ^ = 5. 

C. 

For the speed of light renormalization we present the most detailed results for action 
rather than for since the former, for x > 1, is genuinely anisotropic. In Table XIII 
we list separately regular and tadpole diagram contributions, their gauge invariant sum 
{reg + tad) = (Cq^"* no t.i.) and the fully tadpole improved result (Cq with t.i.), all 
for action and at fixed anisotropy x = 4.0. Cq^^ is independent of ag within numerical 
integration errors which are the most severe when using (^) for nonzero but small masses. 
In Figure 3. we plot Cq^^ both with and without tadpole improvement versus m{°\ Table 
XIV summarizes results for several x values with action and in Figure 4. we plot C^^^ 
versus x 

for fixed a,Mj°^ = 1.0. One sees that tadpole improvement has significant effect 
and causes Cg^^ to switch sign for x > 1- Among other things this allows for a smooth x ^ 1 
limit. 

In Table XV we present results for actions and 5^'. In these isotropic actions nontrivial 
Co comes about because our mass-shell condition p = {iMi, 0) distinguishes between spatial 
and temporal directions once Mi > 0. Table XVI summarizes results for action S^. Here 
tadpole improvement does not decrease the magnitude of the one-loop correction, however, 
for a wide range of mass values it is still true that Cq^^ switches sign for x > 1 ^-iid that the 
X ^ 1 limit becomes smoother after tadpole improvement. 



VI. SUMMARY 

We have carried out one-loop perturbative renormalization of the fermion rest mass Mi, 
wave function renormalization Z2 and the speed of light Cq for a range of highly improved 
actions on isotropic and anisotropic lattices. We find that the dependence of the one-loop 
coefficients on the anisotropy parameter x = O's/O't and on the tree-level mass parameter 
agMi'^ is mild, especially after tadpole improvement of the actions. Furthermore, none 
of the coefficients are particularly large. Mi and Cq exhibit smooth behaviour as one ap- 
proaches the massless, large mass, x 1 a^^d large x limits. This also holds for Z2 if more 



18 



physical combinations such as the difference between continuum and lattice wave function 
renormalizations are considered. The next stage in our program would be to extend the 
present calculations to vertex corrections and to matchings between continuum and lattice 
currents and other multi-fermion operators. 
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APPENDIX A: DEFINITIONS AND FEYNMAN RULES 

In this Appendix we summarize definitions for various terms in the lattice actions and 
present Feynman rules for gauge and quark propagators and for vertex functions. 

Covariant Derivatives Acting on Quark Fields 



^{x) = 1- [U,{x)^{x + a^) - Ulix - a,)^{x - a^)] (Al) 



-^{x) = — [U^{x)^>{x + a^) + UHx - a^)^!{x - a^)] - 2 (A2) 



-^{x) [U,{x)U,{x + a^)^{x + 2a,) - U^x - a^)Ul{x - 2a^)*(x - 2a,)] 

- — \U,{x)^{x^a,) - Ul{x - a,)'^{x - a,)\ (A3) 



\jf -^[x] = \ [U,{x)U,{x + a^)*(x + 2a,) + Ul{x - a,)Ul{x - 2a,)^{x - 2a,) 



- 4 — [U,{x)^{x + a,) + Ulix - a,) ^{x - a,)] + 6 ^{x) (A4) 
Field Strength Tensors 

For the unimproved F,i, of the clover action we use 

^n^x) = . 2 2 23 Ua{x)Ui3{x + aa)U_a{x+aa + af3)U_0{x+a/3), (A5) 
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with {{a, f3)}^^ = {(/i, z/), (z/, -/x), (-/i, -z/), (-z/, /i)} for /i 7^ z/ and ?7_^(x + a^) = ^7^(0;). 
The O(a^) improved field strength tensor of the D234 actions is 



-^{U^{x)F^^{x + af,)Ul{x) + Ul{x - a^)Ff,^{x - a^)U^{x - a^)) - (/x ^ z/) 



u 



6 uf^ ui 



(A6) 



The last term ensures that factors of are correctly removed from those contributions 
to and U^F^^U that end up being four link objects rather than six link ones. In a 

one-loop calculation, however, one can set = 1 everywhere in the definition of the field 
strength tensor and this correction term is irrelevant. 

Both the above covariant derivatives and the field strength tensor are dimensionless. Factors 
of 1/at and l/a^ are inserted explicitly where necessary such as in (|]). 

Gauge Propagator 

The isotropic Symanzik improved gauge action has been discussed quite extensively in the 
literature . Here we summarize formulas for the anisotropic generalization. We start from 
the gauge actions Sq or Sq and add to it a gauge fixing term 



9f 



2a 



a 

1 1 

2ag X 



J2 X^dt{atAt) + J2d,{asA, 



(A7) 



(A8) 



with dfj,Afj,{x) = Afj,{x + — A^{x — a^/2). Equation ( [A8| ) expresses Sgf in terms of 

dimensionless gauge fields a^A^. It is convenient to do so, especially since Sq are already 
in dimensionless form with factors of x 1/x properly put in place. If A^{k) is the Fourier 
transform of (a^A^), the quadratic terms in the gauge action become 



0(0)/,// 



2 ^ 



A^{k)M^,{k) A,{-k 



(A9) 



where 



M 



00 



X 

1 

X 

1 

X 



a. 



a, 



1 



a, 



-k{kj kikj c[{j 



(AlO) 
(All) 
(A12) 
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Moj = M. 



jo = X 



and 



— kokj — kokj qoj 



K - 2 sin(|i). 



The q^i, need to be specified only ioi ji ^ u and one has 



1 



k. 



1 + Y^{kt + k^) 

1 + -k' 
12 ^ 



/i ^ 1/ 



(A13) 



(A14) 



(A15) 



(A16) 



We have inverted the 4x4 matrix M^j, using Mathematica keeping g^,^ generaL For both 
gauge actions, Sq and Sq the free gauge propagator has the structure 



D^,{k) = M. 



Ill/ 



OLgk^ikyX + 



with 



(P)2 _ 



fD{kp,qpa,x) . 



(A17) 



(A18) 



The term proportional to the gauge fixing parameter ag has the familiar form 



a. 



kpku 



{k 



2^2 



X 



an 



(A19) 



with the conversion factor ap^ay/a^at mentioned in (pT]). This factor results because we are 
looking at the propagator for dimensionless gauge fields apAp and because we carried out a 
dimensionless Fourier transform. The second term in ( A17 ) is much more complicated. If 
one writes 



fN^{kp, Qpa, x) — ~ fr^ 
X 

fN{kp,qp„,x) = xfN 
fit^kp, Qpij, x) = X kikj 
In {kp, qpa, X) = X hkj Jn 



(A20) 

(A21) 
(A22) 
(A23) 



one can show that //) and all the fj^ are functions only of {xkoY,kj,gpa with no other x 
dependence or odd powers of kp. We have not shown color indices in the above expressions. 
The gluon propagator is diagonal in color. 
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Quark Propagator 

The inverse free quark propagator for Sjj234, is given by 

1 C -> 

atG^ {k) = 270/^0(^0) sin(A;o) + i — ^'^jKj{kj)fAn{kj) + Q{ko,k) 

X i 



(A24) 



with 



C: 



^0 = 1 + ^-^ cos(/i;o) 



^ 3 



3t 



cos{kj) 



(A25) 
(A26) 



and 



^= X 



2(1 + —jsm (y. 



4i 



6 



sin^(A;o) 



2(l + y)sin^(|)-^sin^(A:,) 



+ atm. 



(A27) 



Propagators for the other quark actions can be obtained by setting the appropriate Cj(t) 
equal to zero. Quark propagators are diagonal in color. 



Vertex Functions 

In deriving the one- and two-gluon emission vertices we have used the method described in 
2^. We list again results only for 5^234 • Those for other quark actions follow trivially. The 



general form for a single gluon emission vertex is 

v,ik',k) ^ -i^,w, -wi- Y.^u,w, 



UfJ, 



(A28) 



where /z is the polarization of the emitted gluon, k' the momentum of the outgoing quark 
and k the momentum of the incoming quark. We suppress the color factor T^^ which should 
multiply each of the above terms. Using the variables 



k 



± 



1 



{k'±k)^ 



(A29) 



one has 



Wn 



W 

3 



Co 
X 

X 

1 

X 



c c 

;i + ^) cos(A;+) 1^ cos(2A;+) cosik^) 

C C 

H ) cos(A;^) — — cos(2A;^) cos(A;^ 

3 3 



:i + ^)sin(A:+ 



a 



At 



6 



sin(2A;(j") cos(fcQ ) 



^1 + ^) sin(/c+) - ^sm{2k^) cos{k- 
3 d 



(A30) 
(A31) 
(A32) 
(A33) 
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and 



W"- 



- sin(2A; • ) cos(A;o ) 77 [5 — cos(2A; • ) — cos(2A;q )] 
2 3 

1 _ _ 1 

- sin(2A;Q ) cos(A; • ) 77 [5 — cos(2A; • ) — cos(2A;q )] 

2 3 

- sin(2/Cj^) cos(/c^) - [5 — cos{2kl[) — cos(2A;;^)]. 
2x ^3 



(A34) 
(A35) 
(A36) 



For the clover action the factor |[5 — cos(2A;^) — cos(2A;^ )] in W'^^ should be replaced by 1. 

For the two-gluon emission vertex we do not present the most general result, but restrict 
ourselves to those terms necessary for the tadpole diagram E*"'^. For instance, the (Jp^pF^y 
term does not contribute to the tadpole diagram. We also omit terms that vanish upon 
symmetrizing between the two gluons. If V^^j^^ik' ,k,qi,q2) stands for the emission vertex 
for gluons of momentum and polarization /ij, with kfj_ = k'^ + qi^^ + ^2,/^, one has 



T/(2) 



j7o 

X 
"2 



^1 + ^)sin(/C(f) - ^C3isin(2A;+)cos(^)cos(^; 



cos{2kl) cos(^) cos(^; 



y{2) ^ 

n 2 X 
1 



7j 



;i + ^) sin(A;^ 



^C3sin(2A:;)cos(^)cos(:|^) 



1 + ^) cos(A:;) - \C, cos(2A;+) cos(%^) cos(%^) 



3 ' ' ^ ' 3 
The color factor for these vertex functions is {T"'^T^^)bc- 



(A37) 



(A38) 



APPENDIX B: IR SUBTRACTIONS 



In this Appendix we list the IR subtraction, JF^^^ of equation (0), necessary to control 
numerical integration of IR divergent integrals. A gluon mass, A/a^, is introduced into D^^ 
by replacing the first factor in (|A17| ) by 



^ fBll 



(A;2)2 k^ fc2 + A2 

The lattice wave function renormalization Z2 must reproduce the same IR divergence struc- 
ture as in continuum QCD. For at one-loop the IR divergence is 

ovr 

^[ 2 - K-1)] ln(A2) m>0 (B2) 

We note that by the m = theory we mean one in which the quark mass has been set to 
zero before taking the limit A ^ 0. This is the usual practice in much of the literature on 
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massless lattice perturbation theory. Alternatively one could take the limit am and 
A — > keeping am > A. Since we want to compare with some of the massless literature 
with improved glue actions (e.g. |T^) we adopt the first definition in this article. In our 
massive calculations we do not go to extremely small masses and have not attempted to 
isolate In(am) contributions. 

For our VEGAS integrations it was convenient to separate the d/dpo derivative in (|6l|) into 
two parts 



_d_ 

dpo 



V 

at 



d 
dpo 



VGV„. 



VGVr, 



den 



VGV 



VGV^. 



VGVr 



den 



{VGV, 



den 



,2 ^^^den' 



(B3) 



Corresponding to the two parts with derivatives acting on the numerator or denominator, 
respectively, we introduce two separate subtraction terms and T^^^. These are ob- 

tained by calculating the self energy diagram in continuum Euclidean perturbation theory 
with an appropriate mass mejj and the mass-shell condition p = {imeff,0). The effective 
mass follows from expanding the lattice integrand in (^) about small k and comparing 
with the continuum calculation |13|. It is given in (|6^. After converting to the dimension- 
less integration variables fc^ of ( |6lD one has for the part with the derivative acting on the 
denominator 



f -4x(x^A:g + 6V4)((A:T-&Vfcg) 

^ I (F + A2) ((A;2)2 + 62^^2^2)2 



with 



P (F + A2) ((P)2 + b^x^k^^ 



(B4) 



k' = xX + j:k] 



2asm 



eff- 



(B5) 



The ^-function imposes a cutoff on J-'f^J^ so that it vanishes identically for k"^ > A"^, where A 
is some number < A < tt. The subtraction term can be integrated analytically to give 



^den^ 

_ as_ 
~ 3^ 



m^ 



> 0, A, A) = / ^eff > 0, A, A) 

A + Vb^ + A2 



A2 

-21n(— ) + 2 In 



+ — + 3A2) + 



A2 



+K - 1) 



, ,A2^ 2A2^ „ „^ A(2A2- 



b^ 
-3621 



b^ 



2A{b^ - 6A2 



V62 + A2 - In 



(B6) 



and 



-pdeni 



m,g = 0,A,A) = — 



■1 + 



Otr. 



A2 



(B7) 



The contribution in (P3|) from the derivative acting on the numerator is 
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+ K - i)x 



which leads to 



F""-(m,^ > 0, A, A) = / ^rr(^, ^eff > 0, A, A) 

J k 



^ , , A + VF+A^\ 4A2 4A 

4 In ( : I — + -rrrVb^ + A2 



; 52 ' 52 



In 



'A + VWTf^\ 2A2 

b J ~ 'V 



and 



F""-(m,^ = 0,A,A) 



as 



— 1 

2 + ^ 



ln( 



A2 
A2 



pden _|_ pnum reproduces the IR divergent logarithms of (|B2|) . 
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TABLES 



action 



comments 



parameters 



— Sq + Sdover 



•S^ — <Sq + Sdover 



S^ — Sq + Sdover 



S'-^ — Sq+ >S£)234 



massive and massless 



massive and massless 



massive and massless 



massless 



c§ = 5/3, cf = -1/12, X = 1 
= 3.648, cf = -0.331, x = 1 
X> 1 

= 5/3, cf = -1/12, X = 1 



C-D _ oil _i_ 
<5 —Oq 



massive and massless 



X> 1 



TABLE I. Combinations of gauge and quark actions considered in this article. 
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Action cS^ 





regular 


tadpole 


reg + tad 


t.i. 


(1) 


71 AD 


/r\\ 




ag 


= 1.0 








0.00 


-1.770(3) 


4.298 


2.528(3) 


-3.001 


-0.472(3) 


0.000 


0.01 


-1.718(5) 


4.266 


2.548(5) 


-2.978 


-0.430(5) 


0.037(6) 


0.05 


-1.547(3) 


4.141 


2.594(3) 


-2.891 


-0.297(3) 


0.152(4) 


0.10 


-1.369(3) 


3.991 


2.623(3) 


-2.786 


-0.164(3) 


0.263(4) 


0.50 


-0.463(3) 


3.030 


2.567(3) 


-2.115 


0.451(3) 


0.737(3) 


1.00 


0.092(3) 


2.260 


2.352(3) 


-1.578 


0.774(3) 


0.948(3) 


2.00 


0.530(3) 


1.511 


2.041(3) 


-1.055 


0.986(3) 


1.050(3) 


5.00 


0.743(3) 


1.096 


1.839(3) 


-0.765 


1.074(3) 


1.077(3) 


10.00 


0.752(3) 


1.075 


1.827(3) 


-0.750 


1.077(3) 


1.077(3) 






ag 


= 0.0 








0.00 


-0.472(3) 


3.001 


2.528(3) 


-3.001 


-0.472(3) 


0.000 


0.01 


-0.429(5) 


2.978 


2.549(5) 


-2.978 


-0.429(5) 


0.038(6) 


0.05 


—0.293(3] 


o cm 
2.891 




— 2.891 


—0.293(3 j 


0.1OD(4j 


0.10 


-0.161(3) 


2.786 


2.626(3) 


-2.786 


-0.161(3) 


0.266(4) 


0.50 


0.454(3) 


2.115 


2.569(3) 


-2.115 


0.454(3) 


0.740(3) 


1.00 


0.775(3) 


1.578 


2.353(3) 


-1.578 


0.775(3) 


0.949(3) 


2.00 


0.987(3) 


1.055 


2.042(3) 


-1.055 


0.987(3) 


1.051(3) 


5.00 


1.076(3) 


0.765 


1.841(3) 


-0.765 


1.076(3) 


1.079(3) 


10.00 


1.079(3) 


0.750 


1.829(3) 


-0.750 


1.079(3) 


1.079(3) 



TABLE II. One-loop mass renormalization for action S for different asM\^ values. Results 
are presented for two choices of the gauge parameter Ug = 1 and Og = 0. OsAfi^^osufe same 
as (reg + tad + t.i.). asM-[^j^^ is defined in (|^) and related to Os-^i^jios^b in (f78|). Where errors 
are not indicated explicitly, they are of 0(1) or less in the last digit. 
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Action 5^' 


Action 


Us Ml 




reg + tad 






reg + tad 




ill 


0.00 




1.480 


-0.397 


0.000 


2.213 


-1.287 


0.000 


0.01 




1.505 


-0.358 


0.035 








0.05 




1.566 


—0.242 


0.136 








0.10 




1.612 


-0.131 


0.228 








0.50 




1.690 


0.367 


0.608 








1.00 




1.611 


0.624 


0.770 








2.00 




1.466 


0.806 


0.859 








5.00 




1.376 


0.897 


0.900 








10.00 




1.371 


0.901 


0.901 









TABLE III. One-loop mass renormalization for actions S^' and . Numerical integration 
errors are at the ±0.006 level for a^M^^^ = 0.01 and of 0(4) or less in the last digit otherwise. 
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Action S 







reg + tad 






reg + tad 








X= 1-0 


X = 2.0 


0.00 




2.710 


-0.479 


0.000 


3.240 


-0.881 


0.000 


0.01 




2.729 


-0.436 


0.038 


3.250 


-0.833 


0.039 


0.05 




2.774 


-0.299 


0.157 


3.260 


-0.675 


0.164 


0.10 




2.799 


-0.164 


0.269 


3.251 


-0.512 


0.288 


0.50 




2.717 


0.463 


0.754 


3.060 


0.290 


0.863 


1.00 




2.466 


0.779 


0.955 


2.824 


0.772 


1.178 


2.00 




2.096 


0.961 


1.026 


2.461 


1.155 


1.381 


5.00 




1.842 


1.014 


1.017 


1.978 


1.371 


1.419 


/r\\ 


X = 3.0 


X = 3.6 


0.00 




3.298 


-1.139 


0.000 


3.286 


-1.243 


0.000 


0.01 




3.302 


-1.092 


0.036 


3.289 


-1.195 


0.036 


0.05 




3.312 


-0.918 


0.167 


3.299 


-1.017 


0.167 


0.10 




3.299 


-0.742 


0.293 


3.287 


-0.835 


0.295 


0.50 




3.110 


0.145 


0.896 


3.106 


0.083 


0.905 


1.00 




2.901 


0.695 


1.244 


2.910 


0.656 


1.262 


2.00 




2,589 


1.161 


1.498 


2.620 


1.151 


1.532 


5.00 




2.106 


1.472 


1.582 


2.158 


1.496 


1.636 


/ r\\ 


X = 4.0 


X = 5.0 


0.00 




3.271 


-1.299 


0.000 


3.233 


-1.405 


0.000 


0.01 




3.271 


-1.255 


0.031 


3.233 


-1.360 


0.031 


0.05 




3.285 


-1.070 


0.167 


3.249 


-1.170 


0.168 


0.10 




3.274 


-0.885 


0.296 


3.239 


-0.980 


0.297 


0.50 




3.099 


0.049 


0.910 


3.076 


-0.018 


0.915 


1.00 




2.911 


0.634 


1.270 


2.903 


0.590 


1.283 


2.00 




2.631 


1.143 


1.547 


2.645 


1.124 


1.572 


5.00 




2.184 


1.505 


1.662 


2.228 


1.518 


1.707 




X = 5.3 


X = 6.0 


0.00 




3.222 


-1.429 


0.000 


3.197 


-1.479 


0.000 


0.01 




3.226 


-1.380 


0.035 


3.204 


-1.426 


0.038 


0.05 




3.238 


-1.194 


0.167 


3.213 


-1.241 


0.168 


0.10 




3.228 


-1.003 


0.296 


3.205 


-1.048 


0.296 


0.50 




3.069 


-0.034 


0.915 


3.053 


-0.065 


0.918 


1.00 




2.899 


0.579 


1.285 


2.891 


0.557 


1.290 


2.00 




2.647 


1.119 


1.577 


2.649 


1.108 


1.586 


5.00 




2.237 


1.520 


1.717 


2.254 


1.522 


1.734 



TABLE IV. One-loop mass renormalization for action for several values of the anisotropy 



X = Qs/at- Numerical integration errors are as in Table III. 
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Action S 





reg + tad 






reg + tad 








X= 1-0 


X = 2.0 


0.00 


2.247 


-1.339 


0.000 


2.740 


-1.927 


0.000 


0.01 


2.268 


-1.289 


0.037 


2.744 


-1.878 


0.030 


0.05 


2.329 


-1.121 


0.153 


2.776 


-1.678 


0.157 


0.10 


2.375 


-0.947 


0.265 


2.785 


-1.473 


0.277 


0.50 


2.451 


-0.043 


0.769 


2.728 


-0.396 


0.858 


1.00 


2.342 


0.510 


1.003 


2.623 


0.320 


1.208 


2.00 


2.082 


0.894 


1.075 


2.404 


0.958 


1.453 


5.00 


1.843 


1.012 


1.021 


1.986 


1.349 


1.455 


/r\\ 


X = 3.0 


X = 3.6 


0.00 


2.788 


-2.250 


0.000 


2,777 


-2.370 


0.000 


0.01 


2.794 


-2.195 


0.033 


2.781 


-2.315 


0.032 


0.05 


2.818 


-1.984 


0.159 


2.805 


-2.099 


0.158 


0.10 


2.824 


-1.762 


0.282 


2.811 


-1.872 


0.282 


0.50 


2.759 


-0.601 


0.883 


2.751 


-0.680 


0.888 


1.00 


2.674 


0.178 


1.262 


2.675 


0.120 


1.275 


2.00 


2.506 


0.901 


1.566 


2.528 


0.870 


1.596 


5.00 


2.115 


1.423 


1.639 


2.166 


1.435 


1.701 


/ r\\ 


X = 4.0 


X = 5.0 


0.00 


2.769 


-2.427 


0.000 


2.732 


-2.544 


0.000 


0.01 


2.768 


-2.377 


0.026 


2.734 


-2.490 


0.029 


0.05 


2.792 


-2.152 


0.159 


2.760 


-2.267 


0.156 


0.10 


2.799 


-1.928 


0.278 


2.768 


-2.030 


0.282 


0.50 


2.742 


-0.722 


0.886 


2.719 


-0.798 


0.891 


1.00 


2.672 


0.090 


1.279 


2.660 


0.033 


1.288 


2.00 


2.536 


0.853 


1.609 


2.542 


0.817 


1.629 


5.00 


2.190 


1.436 


1.729 


2.232 


1.436 


1.779 




X = 5.3 


X = 6.0 


0.00 


2.720 


-2.572 


0.000 


2.701 


-2.620 


0.000 


0.01 


2.724 


-2.516 


0.031 


2.708 


-2.561 


0.033 


0.05 


2.759 


-2.283 


0.166 


2.732 


-2.337 


0.158 


0.10 


2.765 


-2.047 


0.291 


2.740 


-2.099 


0.283 


0.50 


2.717 


-0.810 


0.899 


2.699 


-0.848 


0.894 


1.00 


2.658 


0.022 


1.293 


2.646 


-0.007 


1.291 


2.00 


2.543 


0.809 


1.633 


2.541 


0.791 


1.638 


5.00 


2.242 


1.436 


1.791 


2.256 


1.432 


1.808 



TABLE V. One-loop mass renormalization for action S for several values of the anisotropy 



X = as/at- Numerical integration errors are as in Table III. 
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Action cS^ 





Ug = 1.0 


Ug = 0.0 




( 1 1 

L no t.i. 




f 1 1 

Z2 with t.i. 


f 1 1 

L. no t.i. 

2,a'P0 


^2 jpo ^ith t.i. 


Z2 with t.i. 


0.00 


0.687(3) 


—0.063(3) 


—0.063(3) 


1.195(3) 


0.445(3) 


0.445(3) 


O.Oi 


-2.688(6) 


-3.438(6) 


-3.475(7) 


-2.169(6) 


-2.919(6) 


-2.956(7) 


0.05 


-1.569(5) 


-2.319(5) 


-2.471(6) 


-1.052(5) 


-1.802(5) 


-1.954(6) 


0.10 


-1.041(5) 


-1.791(5) 


-2.054(6) 


-0.524(5) 


-1.274(5) 


-1.537(6) 


0.50 


0.393(3) 


-0.357(3) 


-1.094(4) 


0.905(3) 


0.155(3) 


-0.582(4) 


1.00 


1.112(3) 


0.362(3) 


-0.586(4) 


1.623(3) 


0.873(3) 


-0.075(4) 


2.00 


1.795(3) 


1.045(3) 


-0.005(4) 


2.304(3) 


1.554(3) 


0.504(4) 


5.00 


2.223(20) 


1.473(20) 


0.396(20) 


2.719(20) 


1.969(20) 


0.892(20) 



TABLE VI. One-loop wave function renormahzation for action for different agM^^^ values. 
Results are presented for two choices for the gauge parameter ag = 1 and ag = 0. ^2^'* includes 



both Z, 



(1) 



2 dpo ^^'^ -^2^1/1 defined in equation (82). 



Action ag = 1.0 







^2.dif f with t.i. 


0.00 


0.634(3) 


-0.116(3) 


0.01 


0.668(6) 


-0.082(6) 


0.05 


0.763(5) 


0.013(5) 


0.10 


0.849(5) 


0.099(5) 


0.50 


1.259(3) 


0.509(3) 


1.00 


1.536(3) 


0.786(3) 


2.00 


1.778(3) 


1.028(3) 


5.00 


1.623(20) 


0.873(20) 



TABLE VIL Z2 jpij with In(om) contributions subtracted out for action S in Feynman gauge. 
^i^diff i^ defined in equation (|84|) . 
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Action S^' 







"9 = 


1.0 








"9 = 


0.0 






/r\\ 




7(1) 


with t.i. 


^2 


with t.i. 


^2 jpo ^0 


7(1) 


with t.i. 


^2 


with t.i. 


0.00 


0.239 




-0.230 




-0.230 


0.747 




0.278 




0.278 


0.01 


-3.137 




-3.606 




-3.641 


-2.628 




-3.097 




-3.132 


0.05 


-2.024 




-2.493 




-2.629 


-1.515 




-1.984 




-2.120 


0.10 


-1.502 




-1.971 




-2.199 


-0.993 




-1.462 




-1.690 


U.OU 


n 1 on 




-u.oyu 




— i.iyo 


U.ooy 




n HQ 1 
-U.Uoi 




-U.ooy 


1.00 


0.546 




0.076 




-0.694 


1.054 




0.585 




-0.185 


2.00 


1.165 




0.695 




-0.164 


1.674 




1.204 




0.345 


5.00 


1.568 




1.099 




0.199 


2.079 




1.610 




0.710 












Action S 


C 














«9 = 


1.0 








«9 = 


0.0 










7(1) 


with t.i. 


^(1 
^2 


' with t.i. 


41. no U. 


7(1) 


with t.i. 


7(1) 
^2 


with t.i. 


0.00 


0.145 




-0.355 




-0.355 


0.653 




0.153 




0.153 



TABLE VIII. One-loop wave function renormalization for actions S^' and . Numerical 



integration errors are at the ±0.02 level for a^iVfj = 5.0, at the ±0.006 level for a^Mj = 0.01, 
0.05 and 0.10 and of 0(4) in the last digit for other masses. 
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Action a„ = 1.0 







^24p^ with t.i. 


with t.i. 




^2X0 with 


Z2''' with t.i. 




x = i-o 


X = 2.0 


0.00 


0.789 


-0.023 


-0.023 


0.323 


0.110 


0.110 


0.01 


-2.586 


-3.398 


-3.431 


-3.050 


-3.263 


-3.293 


0.05 


-1.469 


-2.281 


-2.437 


-1.930 


-2.143 


-2.301 


0.10 


-0.942 


-1.755 


-2.024 


-1.406 


-1.619 


-1.887 


0.50 


0.467 


-0.345 


-1.099 


-0.111 


-0.324 


-0.974 


1.00 


1.150 


0.337 


-0.618 


0.398 


0.185 


-0.569 


2.00 


1.766 


0.953 


-0.073 


0.822 


0.609 


-0.147 


5.00 


2.124 


1.312 


0.295 


1.237 


1.025 


0.312 




X = 3.0 


X = 3.6 


0.00 


0.265 


0.172 


0.172 


0.253 


0.190 


0.190 


0.01 


-3.107 


-3.200 


-3.231 


-3.118 


-3.182 


-3.213 


0.05 


-1.988 


-2.080 


-2.240 


-1.999 


-2.062 


-2.222 


0.10 


-1.465 


-1.558 


-1.828 


-1.476 


-1.540 


-1.810 


0.50 


-0.198 


-0.290 


-0.923 


-0.216 


-0.279 


-0.907 


1.00 


0.259 


0.167 


-0.534 


0.228 


0.165 


-0.522 


2.00 


0.594 


0.501 


-0.150 


0.536 


0.473 


-0.147 


5.00 


0.911 


0.819 


0.272 


0.807 


0.744 


0.256 


/ r\\ 


X = 4.0 


X = 5.0 


0.00 


0.249 


0.198 


0.198 


0.242 


0.210 


0.210 


0.01 


-3.123 


-3.174 


-3.205 


-3.130 


-3.162 


-3.193 


0.05 


-2.003 


-2.054 


-2.214 


-2.010 


-2.042 


-2.202 


0.10 


-1.481 


-1.532 


-1.803 


-1.487 


-1.519 


-1.790 


0.50 


-0.223 


-0.275 


-0.901 


-0.233 


-0.266 


-0.889 


1.00 


0.215 


0.164 


-0.517 


0.196 


0.164 


-0.508 


2.00 


0.511 


0.460 


-0.146 


0.473 


0.441 


-0.143 


5.00 


0.758 


0.707 


0.249 


0.678 


0.645 


0.236 




X = 5.3 


X = 6.0 


0.00 


0.241 


0.212 


0.212 


0.238 


0.216 


0.216 


0.01 


-3.131 


-3.160 


-3.191 


-3.133 


-3.156 


-3.186 


0.05 


-2.011 


-2.040 


-2.200 


-2.013 


-2.036 


-2.196 


0.10 


-1.488 


-1.517 


-1.788 


-1.490 


-1.513 


-1.783 


0.50 


-0.236 


-0.264 


-0.887 


-0.239 


-0.261 


-0.882 


1.00 


0.193 


0.164 


-0.506 


0.186 


0.164 


-0.502 


2.00 


0.466 


0.437 


-0.142 


0.453 


0.431 


-0.140 


5.00 


0.661 


0.633 


0.234 


0.631 


0.609 


0.229 



TABLE IX. One-loop wave function renormalization for action S in Feynman gauge. Numer- 
ical integration errors are at the ±0.02 level for agM^^^ = 5.0 and % < 3, at the ±0.006 level for 
a^Mf ^ = 0.01, 0.05 and 0.10 and of 0(4) in the last digit for all other cases. 
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Action a„ = 0.0 





4]L 


^2,L wit^ 


^2^'' with t.i. 




^Xdp^ with t.i. 


with t.i. 




x = i-o 


X = 2.0 


0.00 
0.01 
0.05 
0.10 
0.50 
1.00 
2.00 
5.00 


1.297 
-2.077 
-0.960 
-0.433 
0.975 
1.658 
2.274 
2.635 


0.485 
-2.889 
-1.772 
-1.246 
0.163 
0.845 
1.461 
1.823 


0.485 
-2.930 
-1.931 
-1.518 
-0.593 
-0.111 
0.434 
0.804 


0.833 
-2.539 
-1.419 
-0.895 
0.400 
0.908 
1.332 
1.748 


0.620 
-2.751 
-1.632 
-1.108 
0.187 
0.696 
1.119 
1.535 


0.620 
-2.793 
-1.794 
-1.379 
-0.465 
-0.058 
0.363 
0.822 




X = 3.0 


X = 3.6 


0.00 
0.01 
0.05 
0.10 
0.50 
1.00 
2.00 
5.00 


0.763 
-2.608 
-1.489 
-0.966 
0.301 
0.757 
1.092 
1.408 


0.670 
-2.701 
-1.581 
-1.059 
0.208 
0.665 
0.999 
1.316 


0.670 
-2.742 
-1.745 
-1.332 
-0.426 
-0.037 
0.348 
0.770 


0.747 
-2.624 
-1.505 
-0.982 
0.279 
0.722 
1.029 
1.301 


0.683 
-2.688 
-1.568 
-1.045 
0.215 
0.658 
0.966 
1.237 


0.683 
-2.728 
-1.732 
-1.319 
-0.415 
-0.029 
0.346 
0.750 




X = 4.0 


X = 5.0 


0.00 
0.01 
0.05 
0.10 
0.50 
1.00 
2.00 
5.00 


0.740 
-2.631 
-1.511 
-0.988 
0.269 
0.706 
1.002 
1.249 


0.689 
-2.682 
-1.562 
-1.040 
0.217 
0.655 
0.951 
1.198 


0.689 
-2.722 
-1.726 
-1.313 
-0.410 
-0.026 
0.345 
0.740 


0.729 

-2.642 
-1.522 
-0.999 
0.255 
0.684 
0.960 
1.165 


0.697 
-2.674 
-1.554 
-1.031 
0.222 
0.651 
0.928 
1.133 


0.697 
-2.713 
-1.718 
-1.305 
-0.403 
-0.021 
0.345 
0.724 




X = 5.3 


X = 6.0 


0.00 
0.01 
0.05 
0.10 
0.50 
1.00 
2.00 
5.00 


0.727 
-2.644 
-1.524 
-1.001 
0.252 
0.679 
0.952 
1.148 


0.698 
-2.673 
-1.553 
-1.030 
0.223 
0.650 
0.923 
1.119 


0.698 
-2.712 
-1.716 
-1.303 
-0.401 
-0.020 
0.344 
0.720 


0.723 
-2.648 
-1.528 
-1.005 
0.247 
0.672 
0.937 
1.116 


0.701 
-2.670 
-1.550 
-1.027 
0.225 
0.649 
0.915 
1.093 


0.701 
-2.709 
-1.714 
-1.301 
-0.398 
-0.017 
0.345 
0.714 



TABLE X. One-loop wave function renormahzation for action in Landau gauge. Numerical 



integration errors are as in Table IX. 
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Action ag = 1.0 







^24p^ with t.i. 


with t.i. 


^24po no t.z. 


^2X0 with 


Z2''' with t.i. 




x = i-o 


X = 2.0 


0.00 


0.619 


-0.194 


-0.194 


0.102 


-0.111 


-0.111 


0.01 


-2.755 


-3.568 


-3.597 


-3.270 


-3.482 


-3.511 


0.05 


-1.636 


-2.448 


-2.601 


-2.146 


-2.358 


-2.510 


0.10 


-1.107 


-1.920 


-2.185 


-1.617 


-1.830 


-2.088 


0.50 


0.320 


-0.493 


-1.262 


-0.287 


-0.500 


-1.146 


1.00 


1.034 


0.221 


-0.781 


0.258 


0.045 


-0.728 


2.00 


1.710 


0.898 


-0.178 


0.735 


0.523 


-0.272 


5.00 


2.138 


1.325 


0.304 


1.217 


1.005 


0.274 


/r\\ 


X = 3.0 


X = 3.6 


0.00 


0.038 


-0.054 


-0.054 


0.027 


-0.036 


-0.036 


0.01 


-3.332 


-3.425 


-3.454 


-3.344 


-3.407 


-3.436 


0.05 


-2.209 


-2.302 


-2.452 


-2.221 


-2.284 


-2.435 


0.10 


-1.682 


-1.775 


-2.033 


-1.694 


-1.757 


-2.015 


0.50 


-0.383 


-0.476 


-1.098 


-0.403 


-0.467 


-1.082 


1.00 


0.112 


0.019 


-0.692 


0.078 


0.015 


-0.679 


2.00 


0.499 


0.406 


-0.274 


0.438 


0.375 


-0.271 


5.00 


0.878 


0.786 


0.220 


0.770 


0.707 


0.200 


/ r\\ 


X = 4.0 


X = 5.0 


0.00 


0.023 


-0.028 


-0.028 


0.016 


-0.016 


-0.016 


0.01 


-3.348 


-3.399 


-3.428 


-3.354 


-3.387 


-3.415 


0.05 


-2.225 


-2.276 


-2.426 


-2.231 


-2.264 


-2.414 


0.10 


-1.698 


-1.749 


-2.007 


-1.705 


-1.737 


-1.995 


0.50 


-0.412 


-0.463 


-1.075 


-0.424 


-0.456 


-1.063 


1.00 


0.064 


0.013 


-0.674 


0.043 


0.011 


-0.663 


2.00 


0.412 


0.361 


-0.269 


0.372 


0.340 


-0.265 


5.00 


0.720 


0.668 


0.191 


0.635 


0.603 


0.176 




X = 5.3 


X = 6.0 


0.00 


0.015 


-0.014 


-0.014 


0.013 


-0.009 


-0.009 


0.01 


-3.355 


-3.384 


-3.413 


-3.358 


-3.380 


-3.409 


0.05 


-2.233 


-2.261 


-2.411 


-2.235 


-2.257 


-2.407 


0.10 


-1.706 


-1.735 


-1.992 


-1.709 


-1.731 


-1.988 


0.50 


-0.426 


-0.455 


-1.061 


-0.430 


-0.452 


-1.057 


1.00 


0.039 


0.010 


-0.661 


0.031 


0.009 


-0.658 


2.00 


0.364 


0.335 


-0.264 


0.350 


0.328 


-0.262 


5.00 


0.618 


0.589 


0.173 


0.586 


0.563 


0.168 



TABLE XL One-loop wave function renormalization for action S in Feynman gauge. Numer- 
ical integration errors are as in Table IX. 
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Action ctg = 0.0 





4]L 


^2,L wit^ 


^2^'' with t.i. 




^Xdp^ with t.i. 


with t.i. 




x = i-o 


X = 2.0 


0.00 
0.01 
0.05 
0.10 
0.50 
1.00 
2.00 
5.00 


1.127 
-2.246 
-1.127 
-0.598 
0.828 
1.542 
2.219 
2.646 


0.314 
-3.058 
-1.939 
-1.410 
0.016 
0.730 
1.407 
1.834 


0.314 
-3.097 
-2.095 
-1.679 
-0.756 
-0.274 
0.330 
0.811 


0.612 
-2.758 
-1.635 
-1.106 
0.223 
0.768 
1.246 
1.727 


0.399 
-2.971 
-1.847 
-1.319 
0.011 
0.555 
1.033 
1.514 


0.399 
-3.007 
-2.003 
-1.580 
-0.637 
-0.219 
0.238 
0.784 




X = 3.0 


X = 3.6 


0.00 
0.01 
0.05 
0.10 
0.50 
1.00 
2.00 
5.00 


0.537 
-2.833 
-1.710 
-1.183 
0.116 
0.610 
0.996 
1.376 


0.444 
-2.926 
-1.803 
-1.276 
0.023 
0.518 
0.904 
1.283 


0.444 
-2.960 
-1.959 
-1.537 
-0.600 
-0.194 
0.224 
0.717 


0.521 
-2.849 
-1.727 
-1.200 
0.091 
0.572 
0.932 
1.263 


0.457 
-2.913 
-1.790 
-1.263 
0.027 
0.509 
0.868 
1.200 


0.457 
-2.947 
-1.946 
-1.524 
-0.590 
-0.186 
0.222 
0.693 




X = 4.0 


X = 5.0 


0.00 
0.01 
0.05 
0.10 
0.50 
1.00 
2.00 
5.00 


0.514 
-2.856 
-1.733 
-1.207 
0.080 
0.556 
0.903 
1.210 


0.463 
-2.907 
-1.785 
-1.258 
0.029 
0.505 
0.852 
1.159 


0.463 
-2.941 
-1.940 
-1.519 
-0.585 
-0.182 
0.222 
0.682 


0.504 
-2.867 
-1.744 
-1.217 
0.064 
0.531 
0.859 
1.122 


0.471 
-2.899 
-1.776 
-1.250 
0.031 
0.499 
0.827 
1.090 


0.471 
-2.932 
-1.932 
-1.510 
-0.578 
-0.176 
0.223 
0.663 




X = 5.3 


X = 6.0 


0.00 
0.01 
0.05 
0.10 
0.50 
1.00 
2.00 
5.00 


0.502 
-2.868 
-1.746 
-1.219 
0.061 
0.526 
0.851 
1.104 


0.473 
-2.897 
-1.775 
-1.248 
0.032 
0.497 
0.822 
1.075 


0.473 
-2.929 
-1.930 
-1.509 
-0.576 
-0.175 
0.223 
0.660 


0.498 
-2.872 
-1.749 
-1.223 
0.055 
0.517 
0.835 
1.071 


0.476 
-2.895 
-1.772 
-1.246 
0.033 
0.495 
0.813 
1.049 


0.476 
-2.926 
-1.927 
-1.506 
-0.573 
-0.173 
0.224 
0.654 



TABLE XII. One-loop wave function renormalization for action S in Landau gauge. Numer- 
ical integration errors are as in Table IX. 
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Action 5-^ 




= 4.0 






regular 


tadpole 




no t.i. (reg + tad) 


Cq^^ with t.i. 


asM^ 


ag = 1.0 


C\ C\f\ 

u.uu 


-0.618(3) 


1.474 




0.856(3) 


-0.343(3) 


U.Ul 


-0.616(8) 


1.474 




0.858(8) 


-0.341(8) 


U.U5 


-0.617(5) 


1.474 




0.857(5) 


-0.342(5) 


U.IU 


-0.617(5) 


1.474 




0.857(5) 


-0.342(5) 


\j.O\J 


-0.608(3) 


1.474 




0.866(3) 


-0.333(3) 


1 on 


-0.587(2) 


1.474 




0.887(2) 


-0.312(2) 




-0.547(2) 


1.474 




0.927(2) 


-0.272(2) 


5 nn 


-0.519(2) 


1.474 




0.955(2) 


-0.244(2) 


in nn 


-0.561(2) 


1.474 




0.913(2) 


-0.286(2) 




ag = 0.0 


0.00 


-0.345(3) 


1.199 




0.854(3) 


-0.345(3) 


U.Ul 


-0.351(8) 


1.199 




0.848(8) 


-0.351(8) 


0.05 


-0.347(5) 


1.199 




0.852(5) 


-0.347(5) 


0.10 


-0.347(5) 


1.199 




0.852(5) 


-0.347(5) 


0.50 


-0.333(3) 


1.199 




0.866(3) 


-0.333(3) 


1.00 


-0.312(2) 


1.199 




0.887(2) 


-0.312(2) 


2.00 


-0.273(2) 


1.199 




0.926(2) 


-0.273(2) 


5.00 


-0.244(2) 


1.199 




0.955(2) 


-0.244(2) 


10.00 


-0.286(2) 


1.199 




0.913(2) 


-0.286(2) 



TABLE XIII. One-loop speed of light renormalization for action for different UgM^ ' values 
at fixed anisotropy x = 4.0. Where errors are not indicated explicitly, they are of 0(1) or less in 
the last digit. 
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Action S 





C^^' no t.i 






with t.i. 


Cq no t.i 






with t.i. 






X = 


1.0 






X = 


2.0 




0.00 


-0.039 






-0.019 


0.668 






-0.210 


0.01 


-0.036 






-0.015 


0.671 






-0.207 


0.05 


-0.038 






-0.018 


0.669 






-0.208 


0.10 


-0.044 






-0.023 


0.669 






-0.209 


0.50 


-0.127 






-0.107 


0.654 






-0.223 


1.00 


-0.265 






-0.245 


0.633 






-0.244 


2.00 


-0.509 






-0.488 


0.582 






-0.295 


5.00 


-0.749 






-0.728 


0.432 






-0.446 






X = 


3.0 






X = 


3.6 




0.00 


0.807 






-0.302 


0.842 






-0.330 


0.01 


0.809 






-0.300 


0.843 






-0.328 


0.05 


0.808 






-0.301 


0.842 






-0.329 


0.10 


0.807 






-0.302 


0.842 






-0.330 


0.50 


0.811 






-0.298 


0.849 






-0.322 


1.00 


0.821 






-0.288 


0.867 






-0.304 


2.00 


0.836 






-0.273 


0.899 






-0.272 


5.00 


0.802 






-0.307 


0.907 






-0.264 






X = 


4.0 






X = 


5.0 




0.00 


0.856 






-0.343 


0.880 






-0.363 


0.01 


0.858 






-0.341 


0.882 






-0.361 


0.05 


0.857 






-0.342 


0.881 






-0.362 


0.10 


0.857 






-0.342 


0.881 






-0.362 


0.50 


0.866 






-0.333 


0.893 






-0.350 


1.00 


0.887 






-0.312 


0.918 






-0.325 


2.00 


0.927 






-0.272 


0.969 






-0.274 


5.00 


0.955 






-0.244 


1.031 






-0.212 






X = 


5.3 






X = 


6.0 




0.00 


0.885 






-0.367 


0.893 






-0.374 


0.01 


0.887 






—0.365 


0.896 






—0.372 


0.05 


0.886 






-0.366 


0.895 






-0.373 


0.10 


0.886 






-0.366 


0.895 






-0.373 


0.50 


0.898 






-0.354 


0.907 






-0.360 


1.00 


0.924 






-0.327 


0.935 






-0.332 


2.00 


0.977 






-0.274 


0.992 






-0.275 


5.00 


1.046 






-0.206 


1.073 






-0.194 



TABLE XIV. One-loop speed of light renormalization for action S^. Numerical integration 
errors are at the ±0.008 level for a^M^ ^ = 0.01, at the ±0.005 level for a^ilff ^ = 0.05 and 0.10 



and at the ±0.003 level or less for other masses. 
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Action 


Action 5^' 




C^^^ no t.i. C^^^ 


with t.i. 


C^^^ no 


t.i. C^'^ 


with t.i. 


0.00 


0.000 


0.000 


0.000 




0.000 


0.01 


0.000 


0.000 


0.000 




0.000 


0.05 


0.000 


0.000 


0.000 




0.000 


0.10 


-0.007 


-0.007 


-0.007 




-0.007 


0.50 


-0.102 


-0.102 


-0.098 




-0.098 


1.00 


-0.258 


-0.258 


-0.243 




-0.243 


2.00 


-0.536 


-0.536 


-0.490 




-0.490 


5.00 


-0.809 


-0.809 


-0.730 




-0.730 



TABLE XV. One-loop speed of light renormalization for actions S"^ and S"^' . Numerical inte- 
gration errors are as in Table XIV. 
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Action S 





C}y^ no t.i 






with t.i. 


(1) 

Cq no t.i 






with t.i. 






X = 


1.0 






X = 


2.0 




0.00 


-0.424 






-0.139 


0.213 






-0.301 


0.01 


-0.421 






-0.137 


0.216 






-0.298 


0.05 


-0.424 






-0.139 


0.213 






-0.301 


0.10 


-0.430 






-0.146 


0.211 






-0.303 


0.50 


-0.525 






-0.241 


0.179 






-0.335 


1.00 


-0.683 






-0.398 


0.129 






-0.385 


2.00 


-0.971 






-0.686 


0.033 






-0.482 


5.00 


-1.275 






-0.991 


-0.180 






-0.694 






X = 


3.0 






X = 


3.6 




0.00 


0.320 






-0.388 


0.345 






-0.415 


0.01 


0.324 






-0.385 


0.348 






-0.412 


0.05 


0.321 






-0.387 


0.346 






-0.414 


0.10 


0.319 






-0.390 


0.344 






-0.415 


0.50 


0.305 






-0.404 


0.334 






-0.426 


1.00 


0.286 






-0.422 


0.324 






-0.436 


2.00 


0.256 






-0.452 


0.311 






-0.448 


5.00 


0.167 






-0.541 


0.266 






-0.494 






X = 


4.0 






X = 


5.0 




0.00 


0.356 






-0.427 


0.372 






-0.446 


0.01 


0.358 






-0.425 


0.375 






-0.443 


0.05 


0.356 






-0.426 


0.373 






-0.444 


0.10 


0.355 






-0.427 


0.372 






-0.446 


0.50 


0.347 






-0.436 


0.367 






-0.451 


1.00 


0.340 






-0.443 


0.365 






-0.453 


2.00 


0.335 






-0.447 


0.372 






-0.446 


5.00 


0.311 






-0.471 


0.382 






-0.436 






X = 


5.3 






X = 


6.0 




0.00 


0.375 






-0.450 


0.381 






-0.457 


0.01 


0.378 






—0.446 


0.384 






—0.453 


0.05 


0.376 






-0.448 


0.383 






-0.455 


0.10 


0.376 






-0.449 


0.382 






-0.455 


0.50 


0.371 






-0.454 


0.378 






-0.459 


1.00 


0.370 






-0.455 


0.379 






-0.459 


2.00 


0.379 






-0.446 


0.392 






-0.446 


5.00 


0.396 






-0.429 


0.421 






-0.416 



TABLE XVI. One-loop speed of light renormalization for action S . Numerical integration 
errors are as in Table XIV. 
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FIG. 1. o-s^i^sub '^srsus UgM^^^ both with and without tadpole improvement for action S^. 
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FIG. 2. o-s^i^lub versus x both with and without tadpole improvement for action 
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FIG. 3. C^^^ versus asM[^^ both with and without tadpole improvement for action . 
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Cq versus % both with and without tadpole improvement for action . 
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